Artificial graphene as a tunable Dirac material 
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Artificial honeycomb lattices offer a tunable platform to study massless Dirac quasiparticles and 
their topological and correlated phases. Here we review recent progress in the design and fabrica- 
tion of such synthetic structures focusing on nanopatterning of two-dimensional electron gases in 
semiconductors, molecule-by-molecule assembly by scanning probe methods, and optical trapping 
of ultracold atoms in crystals of light. We also discuss photonic crystals with Dirac cone dispersion 
and topologically protected edge states. We emphasize how the interplay between single-particle 
band structure engineering and cooperative effects leads to spectacular manifestations in tunneling 
and optical spectroscopies. 
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INTRODUCTION 

Graphene is boasting a profound impact in condensed- 
matter science [1-5] and technology [6-8]. It is an unusu- 
ally perfect realization of a 2D semimetal displaying, in a 
wide range of energies, Unearly dispersing conduction and 
valence bands, which touch at the so-called Dirac point. 
Charge neutrahty pins the Fermi energy at the apex of 
the Dirac cones, which are particularly intriguing since 
they imply quasiparticles that behave Hke relativistic el- 
ementary particles with zero rest mass [9]. Dirac cones 
characterize other materials. For instance, they describe 
the chiral conducting surface states that emerge on the 
surface of 3D topological insulators (TIs) [10, 11]. 

In the case of Dirac materials, the standard analysis 
of fermionic systems, based on concepts hke "effective 
mass" and energy gap, requires an extensive reformula- 
tion. The propagation of Dirac fermions shows unusual 
features [9], traceable to the presence of a new physical 
variable, the sublattice-pseudospin degree of freedom. In- 
teractions between Dirac fermions [12] are akin to those 
studied in quantum electrodynamics (QED) and there- 
fore very different from interactions between Schrodinger 
electrons in ordinary metals and semiconductors. 

Dirac conical singularities can actually emerge in any 
2D lattice. Symmetry arguments, in fact, show that they 
should appear regularly at the corners of the Brillouin 
zone (BZ) in lattices with triangular symmetry. There- 
fore Dirac fermion physics and its technological exploita- 
tion should not be considered restricted to the realm 
of "natural" materials — such as graphene and TIs — but 
could also be explored, in principle, in artificial structures 
displaying triangular symmetry. Advantages of such "ar- 
tificial graphene" (AG) systems are likely to be their tun- 
able properties — including lattice constants, hopping en- 



ergies, and inter-particle interaction strength — and the 
possibility to design and realize artificial defects. Ad- 
ditionally, the creation of such structures is particularly 
attractive for the study of quantum phases driven by cor- 
relation effects, which in natural (single-layer) graphene 
seem to emerge only in ultra-high magnetic fields [12]. 

Recent advances have demonstrated the possibility of 
creating AG in diverse subfields of low-energy physics. 
Current methods of design and synthesis include i) 
nanopatterning of ultra-high-mobility 2D electron gases 
(EGs) [13-18], ii) molecule-by-molecule assembly on 
metal surfaces by scanning probe methods [19], iii) trap- 
ping ultracold fermionic and bosonic atoms in honeycomb 
optical lattices [20-22] , and iv) confining photons in hon- 
eycomb photonic crystals [23-26]. Figure 1 summarizes 
these four different approaches and the relevant tunable 
parameters. Here we review the state of the art of this 
emerging multidisciplinary field and offer perspectives on 
possible future developments. 



DESIGNING AND PROBING DIRAC BANDS IN 
ARTIFICIAL LATTICES 

AG structures may not immediately compete with 
graphene for technological applications. However, they 
offer a playground to observe and comprehend physical 
phenomena related to Dirac energy-momentum disper- 
sion relations in regimes that are difficult to achieve in 
natural graphene. In the following we discuss the most 
relevant AG lattices explored so far in which electrons, 
photons, cold atoms, and ions are the confined particles. 
As we elaborate below, these systems have complemen- 
tary physical properties that enable the investigation of 
a wide range of phenomena. 




FIG. 1: Artificial grapiiene structures experimentally obtained by different methods, (a) Scanning electron micro- 
graph of the surface of a nanopatterned GaAs heterostructure. This is achieved by defining first an array of Nickel disks with 
the desired geometry by e-beam nanolithography and then by etching away the material outside the disks by inductive-coupled 
reactive ion shallow etching and by Nickel removal. The bottom panel lists the main tunable parameters of each approach: Vq 
is the lattice depth, d is the lattice constant, and A'' the number of sites. The parameter U {V) is the strength of the on-site 
(nearest-neighbor interaction) repulsion while t is the hopping energy scale. Finally, Tp is the Fermi temperature. Panels (b) 
and (c) refer to molecular graphene systems and optical lattices for cold atoms, respectively. In panel (b) red (black) spheres 
represent the Oxygen (Carbon) atoms of CO molecules while the yellow/orange surfaces represent the electron density in a 
honeycomb pattern. Panel (d) describes photonic honeycomb crystals induced by optical induction methods. An is the change 
of refractive index induced by laser irradiation. The parameter cq is the coupling constant between waveguides, which plays 
the role of the hopping parameter t. There is no Hubbard gap or band gap in this system (only the one that is induced by the 
strain in the form of Landau levels). Panel (c) courtesy of L. Fallani. 



Confining electrons 

In 1970 Esaki and Tsu realized the possibility to en- 
gineer energy bands by artificially modulating the po- 
tential in one direction [27]. They came up with the 
idea of using a semiconductor superlattice and predicted 
the onset of negative differential conductivity. This pio- 
neering work stimulated a large effort worldwide focused 
on band-gap engineering in semiconductor heterostruc- 
tures that, thanks to the refinement of nanofabrication 
techniques, eventually led to the development of lateral 
superlattices — semiconductor systems characterized by a 
2D periodic potential modulation [28]. 

The advent of these artificial crystals with a tun- 
able band-structure has greatly influenced the field of 
2DEGs in modulation-doped semiconductor heterostruc- 
tures. Works in this area started in the late eighties and 
enabled the observation of Weiss oscillations [29], novel 
conductance resonances due to quantization of the elec- 
tron orbits in the 2D pattern in a magnetic field [30], 
chaotic dynamics, and, more recently, led to studies of 
Hofstadter butterfly phenomena [31, 32]. Recently, simi- 
lar effects have been observed in natural graphene where 
the 2D periodic potential (with periodicity of the order 



of ^ 10 nm) was induced by placing it on h-BN [33-35]. 

The ever-increasing toolbox of nano-fabrication meth- 
ods allows today a large flexibility in realizing high qual- 
ity 2D patterns with nanoscale dimensions in semicon- 
ductor quantum structures hosting ultra-high mobility 
electrons. An external potential landscape with honey- 
comb geometry that acts as a lattice of potential wells 
(such as quantum dots) to trap electrons and/or holes can 
be obtained by a combination of e-beam nano lithogra- 
phy, reactive ion etching and deposition of metallic gates. 
The spatial resolution of these techniques can reach val- 
ues of a few tens of nanometer or even below. Further 
improvements in spatial resolution can be obtained by 
bottom-up nanofabrication methods, e.g. by designing 
semiconductor lattices by nanocrystal self assembly [36]. 
The possibility of independently controlling the electron 
density and inter-site distances allows one to tune the in- 
terplay between on-site (U) and nearest-neighbor (V) re- 
pulsive interactions and single-particle hopping (t) , open- 
ing the way to the observation of collective phenomena 
and quantum phase transitions in such AG solid-state 
systems. 

Recent experimental results obtained in honeycomb 
patterns deflned on 2DEGs in GaAs quantum het- 



erostructures offer exciting evidence that AG in semi- 
conductors can be realized in the laboratory [14-17]. 
Available experimental results in the regime U /t ^ 1 
reveal [16] unique low-lying collective excitations, such 
as "anomalous spin waves" , in spectra of inelastic light 
scattering (see also below and Fig. 5a-c). Theoretical 
studies [13, 14, 17] indicate that Dirac bands can be de- 
signed to occur under realistic conditions. For example, 
elementary tight-binding calculations indicate that Dirac 
cones extending for 1 meV can be obtained by tuning 
the quantum dot spacing to 20 nm, a value reachable by 
state-of-the-art top-down nanofabrication methods. 

A completely different route to realize AG with solid- 
state materials has been recently followed in Ref. [19] 
(Fig. 2). These authors succeeded in making AG struc- 
tures with a lattice constant of a few nm by placing 
CO molecules on top of a Cu substrate with the aid of 
the tip of a scanning tunneling microscope (STM) [19] 
(Fig. 2a,b) . By probing the density of states of the con- 
fined electrons via STM measurements (Fig. 2d) and their 
evolution as a function of an applied pseudo-magnetic 
field they proved the formation of Dirac bands with 
the characteristic Landau levels of Dirac fermions [5, 9]. 
While the screening exerted by the bulk states under- 
neath the 2DEG on the Cu(lll) surface makes these 
"molecular graphene" structures not ideal candidates for 
exploring many-body effects, the large versatility in the 
atomic design allows unprecedented local control to em- 
bed, map, and tune the symmetries underlying the 2D 
Dirac equation. In this system, the authors estimate 
U/t ^ 0.5 using known material parameters [19, 37]. 
Within the Cu system, there is room to increase the 
strength of effective interactions by reducing t with larger 
lattices, or else the atomic manipulation scheme may be 
extended to other substrates with lower screening effects 
in a quest to realize other interacting phases (see below) . 

The band structure of molecular graphene can be un- 
derstood by assuming that the superlattice potential cre- 
ated by the CO molecules acts as a weak perturbation on 
the parabolic band that describes the Cu surface state. 
The superlattice potential is most effective at changing 
the parabolic dispersion at the edges of the superlattice 
BZ. Results of such a perturbative calculation [13] are 
shown in Fig. 2c. The superlattice potential hybridizes 
the six unperturbed Cu surface states which lie at the 
corners of the new BZ. The effective Hamiltonian at each 
corner of the BZ is given by a 3 x 3 matrix, which gives 
rise to a doublet (blue and green bands in Fig. 2c) and 
a singlet (red band in Fig. 2c). The doubly degenerate 
state leads to an effective Dirac equation. This "nearly 
free" electron scheme can be generalized to include the 
effect of strain [19] and spin-orbit coupling [38]. Similar 
approaches can be applied to describe the appearance 
of Dirac bands in AG in semiconductors [13, 14]. The 
general symmetries of the triangular lattice uniquely de- 
termine these couplings, which have the same form as 



graphene [39, 40]. In particular, spatially patterning the 
hopping via STM atom manipulation allows the genera- 
tion of both gauge (pseudo) electric and magnetic fields 
in molecular graphene (see Fig. 2e,f). Global changes in 
the lattice constant in molecular graphene add a simple 
scalar potential to the Dirac Hamiltonian equivalent to 
an electrical field which changes the chemical potential 
or "doping" [19]. Local changes to the lattice constant 
engineer a strain introducing a vector potential equiva- 
lent to a large perpendicular magnetic field [19, 41, 42] 
here tunable up to 60 Tesla (Fig. 2f ) . Finally, creating an 
alternating bond structure in the form of a Kekule distor- 
tion was shown (Fig. 2e) to attach mass to the formerly 
massless Dirac fermions, akin to the Higgs field [43-45]. 
For future experiments, the availability of semicon- 
ductors (InAs, InSb, etc.) and metals (Ag, Au, etc.) 
with large spin-orbit coupling creates concrete and excit- 
ing possibilities to explore topological phases of artificial 
matter with these approaches [14, 38, 46, 47]. 



Confining photons 

Photonic crystals, an optical analogue of ordinary crys- 
tals, offer an additional route to design energy disper- 
sion relations with characteristic Dirac points [48]. In 
such crystals, the unusual transmission properties near a 
Dirac point [23-25] were predicted and observed experi- 
mentally. 

The state of the art of photonic crystals operating 
in the microwave frequency range is well described in 
Ref. [49] . In this work the crystal is 2D and composed of 
rows of metallic cylinders, which are arranged to form 
a triangular lattice. Electromagnetic waves propagat- 
ing in such a periodic structure, composed of metallic 
cylinders with radius R = 0.25 a, where a is the lattice 
constant, exhibit a dispersion relation with several Dirac 
points. In the vicinity of a Dirac point, the measured 
refiection spectra resemble the STM spectra of graphene 
fiakes [50, 51]. In a subsequent work [52] extremal trans- 
mission through a microwave photonic crystal and the 
observation of edge states close to Dirac points were also 
demonstrated. The authors of Ref. [52] have shown that 
the transmission through this crystal displays a pseudo- 
diffusive [24] 1/L dependence on the thickness L of the 
crystal. In addition, they measured the eigenmode in- 
tensity distributions in a rectangular microwave billiard 
that contains a triangular photonic crystal. Close to the 
Dirac point there appear states at the straight edge of the 
photonic crystal that represent the artificial counterpart 
of the states at a zigzag edge of natural graphene. Op- 
tical analogues of graphene operating in the microwave 
frequency range have been recently used to simulate 
anisotropic honeycomb lattices and to observe topologi- 
cal phase transitions of Dirac points [53] . 

Since 2007 intensive studies of honeycomb lattices con- 
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FIG. 2: Designer Dirac fermions in molecular graphene. (a) Artificial "molecular" graphene [19] is fabricated via 
atom manipulation, and then imaged and locally probed via scanning tunneling microscopy (STM). Carbon monoxide (CO) 
molecules are individually positioned (blue arrow) with the STM tip into a triangular lattice on a Copper substrate, (b) STM 
topograph of molecular graphene during assembly, showing 2D surface-state electrons repelled from the molecules and guided 
into a honeycomb lattice (black lines), (c) Band structure of molecular graphene calculated using the nearly- free electron (NFE) 
model applied to the periodically-perturbed surface electron gas [13]. The relevant parameters have been chosen to match the 
experimental results in Ref. [19]. A Dirac band crossing appears at the K point of the superlattice Brillouin zone. Similar 
band structures can also be obtained in the tight-binding limit [13, 14] and mapped to the NFE model [19]. Progressively more 
exotic variants of graphene have been fabricated using this method [19] : (d) pristine quasi-neutral graphene exhibits emergent 
massless Dirac fermions (d = 19.2 A, t = 90 meV, t' — 16 meV); (e) graphene with a Kekule distortion (fi = 2^2) dresses the 
Dirac fermions with a scalar gauge field creating mass (0.1±0.02 rrie, rrie being the bare electron mass in vacuum); (f) graphene 
with a triaxial strain distortion embeds a vector gauge field condensing a time-reversal-invariant relativistic quantum Hall 
phase (shown here for a large pseudo-magnetic field of 60 Tesla). In the theory panel, images are color representations of the 
strength of the effective Carbon-Carbon bonds (corresponding to tight-binding hopping parameters t), and the curves shown 
are calculated electronic density of states (DOS) from tight-binding (TB) theory. Insets show gapless and gapped Dirac cones 
matching the experimental data. In the experiment panel, images are STM topographs acquired after molecular assembly 
(100 A field of view, T = 4.2 K), and the curves are normalized tunnel conductance spectra obtained from the associated 
nanomaterial. 



structed by the method of optical induction have been 
performed. The honeycomb structure in Ref. [26] was 
induced by the intensity pattern I{x, y) of three interfer- 
ing plane waves, which is translated into a change in the 
refractive index An through the nonlinearity in a pho- 
torefractive crystal. Such a lattice exhibits several Dirac 
points, which have been termed in Ref. [26] diabolical 
points after M.V. Berry and M. Wilkinson. The paraxial 
evolution of the complex amplitude ^ of a probe beam 
propagating in the lattice is governed by a normalized 



Schrodinger-type equation of the form: 
.9* . „, ^0* 
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where Vq controls the relative values of potential depth 
and nonlinearity strength. The wave dynamics in such 
honeycomb photonic lattices has been extensively stud- 
ied [26] offering evidence of the unique phenomenon of 
conical diffraction around the singular diabolical (Dirac) 
points connecting the first and second bands. These find- 
ings have represented the first experimental observation 
of the phenomenon of conical diffraction, predicted by 
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FIG. 3: Landau levels of photons in strained photonic graphene. (a) Dirac cone in spatial spectrum of unstrained 
piiotonic grapiiene, as obtained by solving Eq. (1) with the Ansatz 'i/{x,y,z) = (j){x,y)exp{l3z). As in Fig. 1, co plays the role 
of the hopping parameter t. (b) Spectrum of inhomogeneously strained photonic graphene, according to a strain tensor that 
gives a constant pseudo-magnetic field. The planes are highly degenerate photonic Landau levels [58]; (c) Eigenvalues in the 
unstrained lattice are listed in descending order. Insets show a microscope image of the photonic lattice (top right) and strong 
diffraction of light through the lattice (bottom left); (d) Eigenvalues in the strained lattice listed in descending order, showing 
clear evidence of Landau levels (labeled by "Ln" with n = —2, —1, 0, -1-1, +2, etc). Microscope image of the strained lattice is 
shown (inset top right) , as well as resulting experimental output showing strong optical confinement in Landau level band gaps 
(bottom left). Courtesy of M. Segev and M. Rechtsman. 



W.R. Hamilton in the 19th century, arising here solely 
from a periodic potential. In addition, "honeycomb gap 
solitons" , residing in the gap between the second and the 
third band, were observed, reflecting yet another special 
property of honeycomb photonic lattices. 

In recent years several theoretical papers were pub- 
lished concerning various aspects of the physics of hon- 
eycomb photonic lattices. These include studies of i) PT- 
symmetry {i.e. symmetry under the combined operations 
of parity and time reversal), ii) nonlinear wave dynam- 
ics, iii) the persistence of the Klein effect, and iv) the 
breakdown of conical diffraction due to nonlinear inter- 
actions [54-56] . Perhaps the culmination of these studies 
has been presented in Refs. [57-59]. Ref. [57] describes 
combined theoretical and experimental work on the cre- 
ation and destruction of topological edge states in "op- 
tical graphene" , where, after the application of uniaxial 
strain (compression), two Dirac points merge resulting 
in the formation of a band gap. Effectively, edge states 
are created (destroyed) on the zig-zag ("bearded") edge 
of the structure. Moreover, the authors of Ref. [57] have 
claimed the observation of a novel type of "bearded" edge 
state, which cannot be explained by the standard tight- 
binding theory, while they can be classified as Tamm 
states lacking any surface effect. This is an example that 
highlights how AG structures might provide insights on 



physics beyond that displayed by natural graphene. A 
second complementary work [58] demonstrates the cre- 
ation of synthetic magnetic fields and "photonic Lan- 
dau levels" separated by bandgaps in the spatial spec- 
trum of the structured dielectric lattice, as illustrated 
in Fig. 3. This is the photonic analog of the relativis- 
tic electron Landau levels observed in strained molec- 
ular graphene [19] (Fig. 2f). Finally, a photonic Flo- 
quet TI based on an AG structure of helical waveguides, 
evanescently coupled to one another, has been proposed 
in Ref. [59]. 2D photonic TIs based on optical spin-orbit 
coupling — achieved through the employment of a meta- 
material composed of split-ring resonators — have been 
proposed in Ref. [60]. 



Confining atoms and ions 

Ultracold atoms and ions represent particularly excel- 
lent arenas for the field of quantum simulation. Effi- 
cient quantum simulators of various systems (like Bose- 
or Fermi-Hubbard models, lattice spin models, etc.) have 
either been realized or are within reach in the time span 
of a few years [61]. 

In recent years, a great deal of attention was devoted 
to simulations involving honeycomb optical lattices. The 



main target of these efforts was not to reproduce the 
physics of natural graphene, but rather to seek out novel 
phenomena that are beyond reach in Carbon-based 2D 
honeycomb crystals. Optical lattices are, for example, 
flexible in the sense that their geometry can be changed 
in situ. Cold atoms in optical lattices, moreover, can 
be forced into regimes of parameters or subject to ex- 
ternal fields, which are difficult to achieve in or un- 
accessible to natural graphene. Examples include the 
regime of ultrastrong spin-orbit coupling [62] and non- 
Abelian gauge fields [63, 64] akin to those that appear 
in the Langrangian of quantum chromodynamics. Note 
that with the advent of molecular graphene and Abelian 
gauge fields [19], proposals now also exist for realizing 
non- Abelian gauge fields in solid-state incarnations [65]. 
Last, but not least, ultracold gases offer novel means to 
control the nature, strength, and range of inter-particle 
interactions [61]. 

The authors of Ref. [21] have pioneered attempts to 
create flexible honeycomb lattices. In this work the first 
realization of an ultracold (^^Rb) Bose gas in a spin- 
dependent optical lattice with hexagonal symmetry was 
reported. The basic structure of this lattice is illustrated 
in Figs. 4a-e. Three linearly-polarized laser beams inter- 
sect at an angle of 120°, yielding local potential minima 
in an hexagonal geometry (Fig. 4a,b) and a local po- 
larization that alternates between cr^ and a^ when one 
goes from one sublattice to the other. As atoms in a light 
field experience a polarization-dependent ac Stark shift, 
the potential at 0-+ and a^ sites is different for different 
atomic Zeeman substates labeled by mp. The angle be- 
tween the z axis and a homogenous static magnetic field 
B, which defines the quantization axis, is denoted by 
(Fig. 4a). According to the local polarization, the two 
interpenetrating triangular lattices forming the hexago- 
nal lattice are denoted as "(T+" and "cr^" lattices. The 
hexagonal lattice can therefore also be regarded as a tri- 
angular lattice with a basis where the atoms occupy the 
(T+ and a^ sites indicated in Fig. 4c for </) = 0° by green 
and red bullets. In contrast, for (f> = 90° the lattice is 
perfectly hexagonal (and the asymmetry between the two 
sublattices vanishes). 

For (j> = 0° , 180° the effective potential felt by an atom 
can be written as: 



V{x) = Vhox(a;) + mFgFfiBBcsix) , 



(2) 



where the polarization of the light field 'P{x) — where 
V{x) = ±1 for pure a polarizations — is mapped onto 
a pseudo-magnetic field B^six) oc —Vhc^{x)V{x)/^B, 
where gp is the Lande g-factor and ^b the Bohr mag- 
neton. The potential consists of a spin-independent part 
Vhexix) of hexagonal symmetry and a state-dependent 
superlattice emerging from the local pseudo-magnetic 
field Bcs{x). The single-particle energy- momentum dis- 
persion relation of an atom evolves as changes from a 
gapped to a gapless one exhibiting Dirac points (Fig. 4e). 



Note that the maxima of the potential at the center of 
each honeycomb cell prevent direct diagonal tunneling 
through these maxima. The interplay between single- 
particle band structure effects and interactions between 
atoms leads to interesting many-body physics [21], which 
will be discussed in the next Section. 

The authors of Ref. [22] have carried out pioneering 
work on cold (^''K) Fermi gases in honeycomb optical 
lattices. They have created, moved, and merged Dirac 
points in a tunable honeycomb lattice. More recently, the 
same authors have studied double transfer through Dirac 
points in a honeycomb optical lattice [66] . More precisely, 
these authors measured the quasi-momentum distribu- 
tion of the atoms after sequentially passing through two 
Dirac points. They observed a double-peak feature in the 
fraction of atoms transferred to the second band, both 
as a function of the band gap at the Dirac points and 
the quasi-momentum of the trajectory. They could not 
observe coherent Stiickelberg oscillations due to the vari- 
ation of the potential gradient over the atomic cloud size. 
The flexibility of the lattice used in these experiments is 
illustrated in Fig. 4f,g. The laser configuration is based 
on three retro-reflected beams, which by control of the 
detuning S, create a checkerboard, a triangular, a hon- 
eycomb lattice, or a lattice of weakly coupled ID chains. 
These fermionic optical lattice techniques [22] were re- 
cently contrasted [67] with complementary methods for 
assembling electron lattices in molecular graphene [19]. 

More work on cold atoms confined in flexible optical 
lattices will be reviewed in the next Section. 

Trapped ions provide perhaps the best quantum con- 
trolled systems available nowadays [68]. However, their 
application for quantum simulations of AG has been so 
far restricted only to small systems. Experimental stud- 
ies on trapped ion systems focused so far on the sim- 
ulation of the Dirac equation encoded in motional and 
internal degrees of freedom of a trapped ion in 3-1-1 and 
1-1-1 space-time dimensions, and its consequences, such 
as Zitterbewegung [69] and the Klein paradox [70] . There 
is, however, a large ongoing effort toward scalable ion 
systems, employing for instance arrays of microtraps or 
self-assembled lattices, where several hundreds of ions in 
a triangular lattice can be stored [71]. 



STRONGLY CORRELATED PHENOMENA AND 

TOPOLOGICAL PHASES IN HONEYCOMB 

LATTICES 

The most interesting phases of matter, including fer- 
romagnetism and superconductivity, emerge in systems 
with a macroscopically large number of degrees of free- 
dom in the presence of interactions. Here, we discuss how 
some prominent broken-symmetry states can be explored 
in artificial honeycomb lattices. 

The great advantage of these artificial crystals is that 
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FIG. 4: Flexible optical lattices for cold atom gases, (a) Experimental setup for a spin-dependent hexagonal lattice [21]: 
three running laser beams intersect at an angle of 120° with respect to each other. The magnetic field B forms an angle </!> with 
the z azis. (b) Resulting scalar potential in the dipole approximation, (c) The total potential for *^Rb atoms in the Zeeman 
state \F = l,mF — +1) is plotted for (j) — 0°,90°, and (j) = 180°. Insets show the corresponding light polarization on the two 
sublattices. (d) Atomic density distributions in the lowest Bloch band of the potentials plotted in c). (e) Dispersion relation 
of the two lowest-energy bands. The colormap illustrates the curvature of the bands. Note that the results for (j) = 180° can 
be obtained from those of <^ = 0° by exchanging minima and maxima of the total potential and atomic density, as well as red 
and green bullets, (f ) An optical lattice of tunable geometry has been created in Ref. [22] by using three retro-reflected laser 
beams. Beams X and Y interfere and produce a chequerboard pattern, while beam X creates an independent standing wave. 
Their relative position is controlled by the detuning S. (g) A large range of lattices can be realized depending on the intensities 
of the lattice beams. They result from the overlap of chequerboard and square-lattice patterns. White (blue) regions denote 
lower (higher) values of the potential energy. The diagram on the right shows the accessible lattice geometries as a function 
of the lattice depths Vx and Vx associated with beams X and X, respectively. The transition between triangular ("T") and 
dimer ("D") lattices is indicated by a dotted line. Dirac points appear in the honeycomb regime ("H.c"), i.e. to the right 
of the dashed line. The limit Vjf ^ Vx,Vy corresponds to weakly coupled, one-dimensional chains ("ID c"). Panels (a)-(e) 
courtesy of M. Weinberg and K. Sengstock. Panels (f)-(g) courtesy of L. Tarruell, T. Uehlinger, and T. Esslinger. 



interactions among electrons, atoms, ions, and photons 
can be engineered and largely tuned. In artificial semi- 
conductor lattices, for example, one can quench the ki- 
netic energy by applying a magnetic field perpendicular 
to the 2D electron system, thereby emphasizing the role 
of long-range Coulomb interactions [13-16]. In an ultra- 
cold atomic gas system, the strength of interactions be- 
tween atoms can be tuned at will by employing Feshbach 
resonances [61, 72]. Finally, significant photon-photon in- 
teractions can be achieved by employing nonlinear optical 
media, thereby offering the possibility to realize strongly 
interacting fluids of light [73] . 



Hubbard correlations and split bands 

The Hubbard model [74, 75] is the cornerstone of the 
physics of strongly correlated systems. The Hubbard 
Hamiltonian encodes a daunting competition between 
two energy scales [76-78]: the kinetic energy t, which 
measures the overlap between wave functions on neigh- 
boring lattice sites, and the interaction energy [/ > 0, 
which measures the strength of on-site repulsions. The 
single-band Fermi-Hubbard Hamiltonian reads: 



n 



-t 






^E^ 



h,-\ni^^ . 



(3) 



Here cj^ (ci,cr) creates (destroys) a fermion with spin 

(J at site i of the honeycomb lattice and fii^^ — c] ^Ci^a 
is the spin-resolved number operator. In the first term, 
the sum is over all pairs {i,j) of nearest-neighbor sites. 
A straightforward generalization of the Hamiltonian (3) 
that applies to bosons exists and it is usually termed the 
"Bose-Hubbard model" [79, 80]. 

When interactions are neghgible {U/t <C 1), 
the ground state of the honeycomb-lattice Fermi- 
Hubbard Hamiltonian is semimetallic and characterized 
by linearly-dispersive massless-Dirac-fermion conduction 
and valence bands touching at two inequivalent points 
in the BZ. In the non-perturbative regime, U/t ^ 1, 
the ground-state phase diagram of this model has been 
extensively studied by means of quantum Monte Carlo 
(QMC) techniques [12]. Early on it was shown that at 
half filling and for U/t > 5 a semimetal-Mott insula- 
tor transition occurs [81]. Recently, a QMC calculation 
has demonstrated [82] the existence of a gapped anti- 
ferromagnetic phase for U/t > 4.3 and presented evi- 
dence for a gapped spin liquid phase in the range of cou- 
plings 3.5 < U/t < 4.3. The findings of Ref. [82] have 
been recently addressed in Ref. [83], by QMC simula- 
tions of the same model in larger clusters (containing up 
to 2592 sites) , finding very weak evidence of a spin liquid 
phase. 

An "extended" Hubbard model involving an additional 
interaction term describing nearest neighbor repulsions of 



strength V has been studied analytically by employing 
renormalization-group (RG) techniques in the limit of a 
large number N{ of fermion flavors [84, 85] (we remind 
the reader that A^f = 4 for electrons in graphene). It has 
been found [84, 85] that sufficiently large values of V/t 
stabilize charge-density-wave phases over semimetallic or 
Mott insulating phases. Numerical RG calculations [86] 
have qualitatively confirmed these results but also dis- 
covered that second-neighbor repulsions favor topologi- 
cal states with spontaneously broken time-reversal sym- 
metry (i.e. quantum spin Hall phases) over charge- or 
spin-density waves. 

Finally, it has been recently shown [87] that, when the 
Fermi energy is moved away from the Dirac point and the 
system is doped to the vicinity of a van Hove singularity 
in the density of states (such as the one that occurs at the 
M point of the honeycomb-lattice band structure [5, 12]), 
repulsive short-range interactions favor d + id chiral su- 
perconductivity over many other competing orderings. 
An alternative route to non-conventional superconduc- 
tivity in AG and related compounds relies on the sup- 
pression of long-range forces, while keeping a large value 
{U/t ^ 1) of the Hubbard on-site repulsion [88, 89]. 

All these predictions call for future experimental ef- 
forts. To this end we note that strong correlations leave 
deep scars on the excitation spectrum of a many-body 
system, which can be probed by a variety of spectro- 
scopic tools. The simplest example is represented by a 
gapped collective mode between "Hubbard split bands" . 
In the strongly correlated or "atomic" U/t 3> 1 limit the 
Hubbard model (3) displays two bands [75, 90], which 
are split by the energy cost U of having two fermions 
with antiparallel spin on the same site. In the atomic 
limit it is therefore natural to expect a gapped collective 
mode in which particles belonging to the lower Hubbard 
band are cooperatively promoted to the upper Hubbard 
band. Such a mode has been observed in a Bose-Einstein 
condensate in a deep 3D optical lattice [91] and, more 
recently, in an artificial honeycomb lattice realized by 
nanopatterning the surface of a GaAs semiconductor [16] 
(see Fig. 5). This AG system displays also an anoma- 
lous "spin-flip" mode (Fig. 5c) that could arise from the 
removal of the sublattice-pseudospin degeneracy due to 
inter-site Coulomb repulsions [92] . Explorations of Hub- 
bard excitations and anomalous spin-flip modes can be 
extended to molecular and photonic AG once interactions 
in these systems are turned on. 

As we have already stated above, a plethora of many- 
body effects and lattice models can be very effectively 
simulated by employing cold atom gases [61]. For exam- 
ple, in Ref. [21] combined effects of the lattice and inter- 
atomic interactions led to a forced antiferromagnetic Neel 
order, in which two spin-components localize at differ- 
ent lattice sites. Coexistence of Mott-insulator- and 
superfluid-type orders leads to the formation of a forced 
supersolid. Next-nearest-neighbor tunneling seems to 
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FIG. 5: Many-body effects in artificial graphene. (a) Resonant inelastic light scattering spectra (at B = 5.48 Tesla and 
T = 1.7 K) of a high-mobility 2D electron gas in the presence of a honeycomb lateral superlattice [16] showing a weak peak 
in the excitation spectrum due to the ordinary cyclotron mode and a much stronger peak due to the Hubbard mode, (b) 
Evolution of the energies of the cyclotron mode (black filled circles) and of the Hubbard mode (at frequencies a;HB) (red filled 
squares) at T = 1.7 K. The black dashed line is a linear fit to the data, Tiloc ~ TieB/{m*c), with m* — 0.067 rric. The red 
dashed line is a fit of the type Tilohb = a\/-B[Tesla] with a ~ 2 meV. (c) Energies of two spin collective modes: the ordinary 
spin- wave mode (black filled circles) and an anomalous "spin-flip" mode (red filled triangles). The black dashed line is a linear 
fit to the data, Esw ~ IfflA'sBx, with \g\ = 0.42. Representative excitation spectra of both spin modes at two different laser 
energies are reported in the inset (Bt is the total magnetic field and 9 the tilt angle of the sample with respect to the magnetic 
field). The black filled squares label the splitting A between the two spin modes, (d) Superfluid-to-Mott-insulator transition 
in a spin-dependent honeycomb optical lattice [21]. The visibility of the interference fringes in a time-of-flight experiment is 
plotted as a function of the lattice depth (in units of the recoil energy). At the superfiuid-to-Mott- insulator transition the 
visibility drops down considerably. Data labeled by blue symbols refer to the superfluid-to-Mott-insulator transition in a spin- 
independent hexagonal lattice, i.e. rrip = in Eq. (2). The two other curves refer to finite values of mp, as indicated in the 
inset. Panels (a)-(c) courtesy of M. Gibertini (experimental data from Ref. [16]). Panel (d) courtesy of C. Olschlager and P. 
Windpassinger (experimental data from Ref. [21]). 
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play a role in the physics described in Ref. [21], a fact 
that paves the way for the realization of the famous Hal- 
dane model [93] in this system. Later on, the authors 
of Ref. [94] reported the observation of a quantum phase 
transition to a multi-orbital superfluid phase in an opti- 
cal lattice. In this unconventional superfluid, the local 
phase angle of the complex order parameter is contin- 
uously twisted between neighboring lattice sites. The 
nature of this twisted superfluid quantum phase is an 
interaction-induced admixture of the p-orbital contribu- 
tions favored by the graphene-like band structure of the 
hexagonal optical lattice used in the experiment. 

In Ref. [95] various forms of frustrated classical ferro- 
magnetism have been studied by transforming the lat- 
tice geometry from square to triangular to an array of 
linear chains and back. Finally, in a recent experimen- 
tal work [96] novel instances of quantum magnetism of 
ultracold fermions have been reported. In this work 
short-range magnetic order has been achieved by load- 
ing two-component Fermi gases in either a dimerized or 
anisotropic simple cubic optical lattice. Flexible kagome 
lattices have been recently studied in Ref. [97]. 



Long-range interactions 

Electrons in nanopattcrned 2DEGs [13-18] offer a nat- 
ural system to study correlation effects in the presence 
of long-range Coulomb interactions. The applicability 
of the Hubbard model ([74, 75]) to describe these sys- 
tems when the Fermi energy is at the Dirac point [5, 12] 
(the density of states vanishes) is questionable because 
the 1/r long-range tail of the Coulomb interaction is not 
screened. 

In recent years, several efforts have been made [98- 
101] to describe electrons moving in a honeycomb lat- 
tice and interacting through the non-relativistic Coulomb 
force. In this case the strength of interactions is measured 
by the dimensionless parameter [5, 12] Ofcc = e^/i^^vp), 
where e is the absolute value of the electron's charge, e is 
a suitably-defined dielectric constant, and v-p is the Fermi 
velocity. This parameter, which formally resembles the 
QED fine-structure constant a = e^ /{he) (where c is the 
speed of light), plays the role of the U /t coupling constant 
in the Hubbard model [74, 75]. 

A excitonic insulating phase has been predicted to 
occur spontaneously [99-101] at a critical value of the 
Coulomb coupling constant aco — 1.1 for A^f = 4 
fermion flavors. Note that electrons in natural suspended 
graphene are characterized by a^c ~ 2.2 (since, in this 
case, e ^ 1 and v-p ^ 10^ m/s). Weak-field magneto- 
transport experiments have been carried out [102, 103] 
in high quality suspended samples (with mobilities of 
the order of ^ 10^ cm^ V^^s^^), with carrier density 
fluctuations and temperatures as low as 10* cm~^ and 
1 K, respectively. No experimental evidence of a gapped 



phase has been reported so far in nearly-neutral natural 
graphene [102, 103]. The observed behavior of graphene 
is more consistent with the existence of a strong renor- 
malization of the Fermi velocity [102] and an RG flow to- 
wards a weakly coupled phase (see, however, Ref. [104]). 
When one dopes the system away from the neutrality 
point, screening kicks in and the electron fluid in single- 
layer graphene behaves as a Fermi liquid, albeit with a 
number of intriguing twists [12, 105, 106]. 

Electrons in artiflcial honeycomb lattices [13-16] there- 
fore offer the unique opportunity to study strong correla- 
tion effects, which are beyond reach in natural graphene. 
The role of long-range electron-electron interactions in 
artificial graphene has been studied by means of density- 
functional theory in Ref. [107]. Electron-electron interac- 
tions have been demonstrated [107] to shift the threshold 
for the emergence of isolated Dirac points to larger well 
depths than found without interactions [13, 14]. This ef- 
fect is particularly pronounced when the number of elec- 
trons per well is increased. 

Long-range interactions can nowadays also be stud- 
ied in the realm of atomic physics. Recent experimental 
advances in cooling atoms with permanent dipole mo- 
ments and polar molecules have indeed made it possible 
to study quantum gases with dipolar interactions decay- 
ing like 1/r^ at large distances. These systems have at- 
tracted a large fraction of theoretical and experimental 
interest: detailed information can be found in recent re- 
view articles [108, 109]. 



SUMMARY AND PERSPECTIVES 

In summary, in this Article we have reviewed recent 
progress in the creation of artiflcial graphene focusing on 
nanopatterning of ultra-high-mobility two-dimensional 
electron gases in semiconductors, molecule-by-molecule 
assembly via scanning probe methods, conflning photons 
in dielectric crystals, and optical trapping of ultracold 
atoms in crystals of light. Lattices of superconducting 
circuits [110] may also offer further opportunities to study 
strongly correlated phases of light in honeycomb struc- 
tures. Fully tunable plasmonic analogues of graphene 
can be realized in two-dimensional honeycomb lattices of 
metallic nanoparticles [111]. 

The occurrence of fragile interaction-induced broken- 
symmetry states in artiflcial graphene might be pre- 
empted by external unwanted random potentials, which 
induce inhomogeneities in the particle-number distribu- 
tion. Controlled sources of disorder can, however, be ex- 
ploited to investigate the interplay between disorder and 
interactions [112] and the occurrence of "Bose glass" [79] 
phases. 

Finally, we foresee useful applications of artificial 
graphene for the realization of new topological phases of 
matter and the engineering of Abelian and non-Abelian 
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gauge fields. 
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